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ABSTRACT 

We present a detailed analysis of spectral line profiles in Type la supernova (SN la) 
spectra. We focus on the feature at ^3500 - 4000 A, which is commonly thought to 
be caused by blueshifted absorption of Ca H&K. Unlike some other spectral features 
in SN la spectra, this feature often has two overlapping (blue and red) components. 
It is accepted that the red component comes from photospheric calcium. However, it 
has been proposed that the blue component is caused by either high-velocity calcium 
(from either abundance or density enhancements above the photosphere of the SN) or 
Sill A3858. By looking at multiple data sets and model spectra, we conclude that the 
blue component of the Ca H&K feature is caused by Si II A3858 for most SNe la. The 
strength of the Sill A3858 feature varies strongly with the light-curve shape of a SN. 
As a result, the velocity measured from a single-Gaussian fit to the full line profile 
correlates with light-curve shape. The velocity of the Ca H&K component of the profile 
does not correlate with light-curve shape, contrary to previous claims. We detail the 
pitfalls of assuming that the blue component of the Ca H&K feature is caused by 
calcium, with implications for our understanding of SN la progenitors, explosions, 
and cosmology. 

Key words: line: identification - line: profile - supernovae: general - supernovae: 
individual: SN 2010ae - supernovae: individual: SN 2011fe 



1 INTRODUCTION 

The spectral-energy distribution (SED) of a Type la super- 
nova (SN la) near maximum brightness is relatively simi- 
lar to that of a hot star. An SN SED is predominantly a 
black body with line blanketing in the ultraviolet. There are 
also prominent spectral features associated with absorption 
and emission from elements primarily generated in the SN 
explosion. These features typically have broad P-Cygni pro- 
files, although overlapping lines can produce larger and more 
complicated profiles. 

The exact SED of a SN la depend s on the velocity an d 
density structure of the SN ejecta (e.g.. [Branch et al1ll985h . 
Since broad-band filters sample portions of the SED, and 
measurements in such filters are used to determine SN dis- 
tances to ultimate l y measure cosmolo gical parameters (e.g., 
IConlev et al.ll201ll : lsiMuki et al.ll2012l l. understanding SN la 
spectral features is important for precise cosmological mea- 
surements. 

SN spectra also provide detailed information about 
the SN explosion, progenitor compos ition, circumbinary 
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environment, and reddening law (e.g., Hofiic h et al 
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Ropke et alj l2012f ). Furthermore, there is evidence that one 
can estimate the intrinsic colour of SNe la, and thus improve 
distance measurements through a better estimate of the 
dust reddening by m easuring the ejecta velocity of SNe la 
(|Foley fc Kasenlboilh . Ejecta velocity is measured from the 
blueshifted position of spectral features. For both cosmology 
and SN physics, it is important to have a precise understand- 
ing of SN la SEDs. 

At optical wavelengths, the two most prominent fea- 
tures in a maximum-light spectrum of a SN la are at 
~3750 and 6100 A, respectively. The latter is thought to 
be from Sill AA6347, 6371 (g/- weighted rest wavelength of 
6355 A), and is the hallmark spectral feature of a SN la. 
The former, at rest-frame wavelengths of ~3500 - 4000 A, 
is generally attributed to blueshifted absorption from Ca 
H&K (g/- weighted rest wavelength of 3945 A). However, 
the line profile of this feature is complicated, often times 
displaying shoulders, a flat bottom, a "split" profile, and/or 
two distinct absorption components. There is broad con- 
sensus that the red component of the profile is from Ca 
H&K at a "photospheric" velocity, i.e., a velocity simi- 
lar to that of the ejecta at close to the r = 2/3 sur- 
face, which is typically about 12000 kms~^ near maxi- 
mum light. However, there is no clear consensus to the 
origin of the blue component. Previous studies have at- 
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tributed the blue absorption component to either "high- 
velocity" (HV) Ca HfcK absorption (^18,000 kms~^ e.g 
Hatano et al.ll 19991: iGaravini et al.ll2004l: [Branch et al.l 2005 



20071: IStanishev et all l2007l: IChornock fc Fihppenkol I200S 



Tanaka et al.l |2008| . I2OIII : IParrent et all 120121 '). where the 
absorption comes from a region at high velocity within 
the SN ejecta that has high-density calcium, and to Sill 
AA3 854, 3856, 2863 (g.f- weighted r est w avelength of 38 58 A; 
e.g.. iKirshiier et al.lll993l: iHoflichlliggll: iNugent et al Ill997l : 



iLentz et al.l I2OO0I : IWang et al.l l2003l : lAhavilla et al.l |2007[ ). 

Since calcium and silicon produce the strongest features in 
SN la spectra near maximum brightness, both interpreta- 
tions are worth investigation. For convenience, we will gen- 
erally refer to this feature as the "Ca H&K feature." 

There are several cases of clear HV material in 
SNe la. Observations showing mu ltiple components to 



the Sill A6355 lin e profile (e.g., iMazzah et al.l lioOSbl: 



Altavilla ct al.' '2007*; 'Garavini et al 200?': 'St anishev et al I 
2007; Wang ct al. 2009b : Foley et a. 1. 2012^) or strong 
and quickly varyin g HV Ol A7774 (|Ahavilla et al.l l2007l : 
iNugent et al.ll2011^ are perhaps the cleanest way to detect 
HV material since there are no other strong lines just blue- 
ward of Sill A6355 and 01 A7774. Other detections have 



spectropolarimetry (e.g., Hatano et al. 
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Mazzali et al.l l2005al). however there are several subtleties 



to this feature. 

HV features must be caused by abundance and/or den- 
sity enhancements in layers of the ejecta above the SN 
photosphere. Two distinct "layers" of material within a 
smooth density profile (i.e., an abundance enhancement) 
would necessarily be caused by the explosion, and obser- 
vations of HV features could ther efore restri c t the possi- 
ble explosion models. However, Mazzali et al.l (|2005bl ) sug- 
gested that abundance differences alone cannot reproduce 
the strength of the HV features, and therefore there must be 
a density enhancement. Density enhancements may be either 
caused by the explosion causing over-dense blobs or shells 
of material or by sw e eping up circumbinary material (e.g., 
Gerardv et"aLl |2004| : iMazzali et"al1 l2005al : lOuimbv etaU 



20061 ). Spectropolarimetric observations have indicated that 



HV Ca NIR triplet features are pro b ably caused from 
the explosion (e.g., iKasen et al.l l2003l : IWang et al.l l2003l : 
IChornock fc Fihppenkol 20081 ). Because of its wavelength, it 
is difficult to obtain high-quality spectropolari metric mea- 
surem ents of the Ca H&K feature. None the less. I Wang et al.l 
l|2003l ) was able to make such a measurement, and the po- 
larization spectrum suggested that the blue component of 
the Ca HfcK feature was from Sill A3858 for SN 2001el. 

Using the large CfA sample of SN la s pectra 
iBlondin et al.l l2012bl ). iFolev. Sanders, fc Kirshiieil (|201ll ) 
determined that the velocity of Sill A6355, vsui, and the ve- 
locity of the red component of the Ca HfcK feature, ncaH&K, 
at maximum light correlated with intrinsic colour, but did 
not correlate with light-curve shape (and thus luminosity). 
However, they did not find statistically significant evidence 
of a linear correlation between the pseudo-equivalent width 
of the Ca H&K feature and intrinsic colour. Using S DSS-II 
Super nova Survey and Supernova Legacy Survey data. lFtalevI 
(|2012l ) confirmed these trends with high-redshift SNe la. 



They also noted a slight {2A-a significant) trend between the 
ma ximum-light t>caH&:K ( fria, h^k) and h ost-galaxy mass. 

iMaguire et"all l|2012l . hereafter IM12| ) presented a sam- 
ple of maximum-light low-redshift SN la spectra obtained 
with the Hubble Space Telescope (HST). After various qual- 
ity cuts, the sample consisted of 16 spectra of 16 SNe la. 
These spectra covered the Ca HfcK feat ure, but did not 
cov er wave l ength s near Si II A6355. Unlike .Folev et all (|201ll ) 
and iFolevI (|2012l ). which presumed that the red component 
of the Ca HfcK feature was representative of photospheric 
calcium (and thus the wavelength of the m aximu m absorp- 
tion of this component represented tina h^k). IM12| fit a single 
Gaussian to t he en tire profile to measure ncaH&K- Among 
other claims, IM12| reported a linear relationship between 
'^CaH&K and light-curve shape (3.4-(t significant). Further- 
more, they claimed that after correcting for the relation be- 
tween light-curve shape and wcaH&K, there is no correlation 
between iicaH&K and host-gala:xy mass. 

In this paper, we examine the claims of IM12| with par- 
ticular scrutiny to the details of the Ca H&K profile. In 
Section[21 we re-examine the M12 sample. We confirm a dif- 
ference in the Ca HfcK line profile for SNe la with different 
light-curve shapes, but show that the difference is primar- 
ily in the blue component. We also conclude that single- 
Gaussian fits to the Ca HfcK feature give biased, unphysi- 
cal velocity measurements. In Section (3] we provide simple 
models of calcium and silicon features in a SN la spectrum. 
Trends in the spectra indicate that the blue component of 
the Ca H&K feature is likely from Sil l A38 58. In Section H 
we perform further analysis with the sample, finding 

systematic biases in single- Gauss ian velocity measurements 
appear to be present in the IM12| analysis. In Section O we 
re-examine the CfA spectral sample, providing further evi- 
dence that (1) the blue component of the Ca HfcK feature 
is from Sill A3858 absorption and (2) there is no evidence 
for a correlation between vcn. h&k and light-curve shape. Fi- 
nally, we examine the spectra of the well-observed SN 201 Ife 
and SN 2010ae, a very low-velocity SN lax, in Section[6l Al- 
though one cannot uniquely claim that the blue component 
of the Ca HfcK profile is from Sill A3858 for SN 201 Ife, it 
must be the the case for SN 2010ae. We discuss implications 
of this result and conclude in Section [7| 



2 THE CA H&K LINE PROFILE 

As noted above, the spectral feature at rest-frame wave- 
lengths of ~3500 - 4000 A in SN la spectra often has struc- 
ture such as shoulders and multiple components. The main 
absorption is thought to be from Ca HfcK (from the photo- 
sphere and possibly from a HV component) and Sill A3858. 
We will refer to this feature as the Ca HfcK feature, although 
there may be additional species contributing to it. 

In this section, we will examine the Ca H&K feature in 
detail. To do this, we use the lM12l spectra. After testing their 
claims of differences in the profile shape with light-curve 
shape, we examine the different results one gets depending 
on t he m ethod of fitting the line profile. 

IM12| suggested that «caH&K depends on the light-cure 
shape (and therefore pe ak luminosity) of the S N. Using 
the WISERep database jYaron fc Gal- Yamllioi^ ) . we ob- 
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Figure 1. Median spectra from the M12 sample. The black curve 
is the median spectrum from their full sample. In the top panel, 
the blue and red curves represent the median spectra taken from 
early and late subsamples, respectively, while in the bottom panel, 
they represent low and high stretch (corresponding to low and 
high luminosity) subsamples, respectively. 



tained most of th e spectra presented bv lMl j^l . We exclude 
all SNe that IM12| do not use in their final analysis, includ- 
ing PTFlOufj, which only has a redshift determined by SN 
spectral feature matc hing. In total, there are 14 spectra of 
14 SNe la in the final IM12| sa mple. In Figure [T] we present 
median spectra from the lMld sample. The Ca H&K feature 
has two clear minima (at ~3720 and 3800 A, respectively) in 
the median spectrum. These data appear to be an excellent 
sample for studying the Ca H&K profile shape. 

We also generated median spectra for subsets of the full 
sample. First, we split the sample by phase. Since the ve- 
locity of SN features typicaUy decreases monotonically with 
time because of the receding photosphere, one expects lower 
velocity features at later times. The phase-split median spec- 
tra do not appear to be significantly difi'erent from each 
other or the median sp ectrum from the full sample. This is 
likely the resuh of the IM12| sample having a very narrow 
phase range. 

We also split the sample by light-curve shape. W e split 
the sample by s = 1.01 to match what was done bv Im12| . 
Here, we see the same result that IM12| found and shows 
in their Figures 5 and 7. Namely, the low-stretch (corre- 
sponding to faster-declines and lower luminosity) SNe have 
narrower, seemingly lower velocity features than those of 
high-stretch SNe. 

Given the above difference, it is worth a detailed look 
at the line profiles. Despite coming from P-Cygni profiles, 
the line profiles appear to be similar to the sum of two 



^ The spectra of SNe 2011by and 2011fe were not included in the 
database, and are therefore excluded from our analysis. But we 
do not expect including them in our analysis would change our 
results much. 



Gaussians, and performing such a fit resulted in excellent 
matches to the profiles. The absorption component of a P- 
Cygni profile is very similar to a Gaussian, so using Gaus- 
sians to fit the absorption is a reasonable choice. In Fig- 
ure [2] we display the median spectra for the full sample and 
the low/high-stretch subsamples. We also display the best- 
fitting double- Gaussian fits (after removing a linear pseudo- 
continuum) to each line profile. For each case, we performed 
a six-parameter fit, allowing the centroid, width, and height 
of each Gaussian to vary. The centroid of each Gaussian 
corresponds roughly to the characteristic velocity of that 
component. Similarly, the width of each feature corresponds 
to the velocity-width of the absorbing region for that fea- 
ture. Finally, the height of each feature is roughly related to 
the amount of absorbing material at a given velocity. The 
six-parameter double- Gaussian fits to each profile are rep- 
resented by the blue lines in Figure (2] 

We also fit the low/high-stretch subsamples with two 
parameters fixed and four allowed to vary. The centroid 
and width (the parameters related to velocity) of the red- 
der Gaussian was fixed to match the best-fitting values 
for the full-sample median spectrum, and the remaining 
parameters (all parameters for the bluer Gaussian and 
the height of the redder Gaussian) were allowed to vary. 
These fits are represented by the red lines in Figure (2] 
Visually, the six-parameter fit is not a significantly bet- 
ter representation of the data than the four-parameter fit. 
The reduced decreases by 0.10 and 0.06 when changing 
from the six-parameter to the four-parameter fit for the 
low and high-stretch subsamples, respectively. That is, the 
four-parameter fit has a smaller reduced than the six- 
parameter fit (although only marginally smaller), and thus, 
the subsamples and the full sample are completely consistent 
with all h aving the same velocity for the red component. 

IM12| argued that the difference in the red edge of the 
Ca H&K line profile was evidence that the subsamples have 
different ejecta velocities. But we have shown that simply 
varying the height of the redder Gaussian (and the bluer 
Gaussian) are sufficient to produce the red edge of the pro- 
file. That is, the apparent different in the red edge can be 
explained by different line strengths rather than different 
line velocities, and thus a difference in the red edge is not 
sufficient to distinguish different velocity features. 

We also attempted to fix the parameters of the bluer 
Gaussian, but that did not result in good fits. From these 
tests, we see that (1) the red component does not necessarily 
have a different centroid (and thus velocity) for the two sub- 
samples and (2) the blue component does have a different 
centroid. 

We now turn to the difficulty of reducing these pro- 
files to a single parameter, namely velocity. There have been 
two approaches to measure velocities. The first fits a single 
Gaussian to a line profile and ascribes the centroid of the 
Gaussian to the velocity of the fe ature. This method is used 
by many studies, including IM12| . The alternative is to mea- 
sure the wavelength of maximum absorption (usually after 
some smoothing) to represent t he velocity of the fea ture. 
This is the method described bvlBlondin et al.1 (|2006l ) and 
used bv lFolev et all |2oi3) and lFolevI (I2OI2I ). Although there 
are many arguments to use either method, we will focus on 
the potential systematic errors of using these methods when 
a feature has multiple components like the Ca H&K feature. 
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Figure 2. Median spectra from the iMl j sample. Tlie left, center, and right panels show the median spectra from the full sample, the 
low-stretch subsample, and the high-stretch subsample, respectively. The blue and red curves are double- Gaussian fits to the data. The 
blue curve in the left panel and the red curves in the middle and right panels are the best fits with all parameters allowed to vary. 
The blue curves in the middle and right panels represent fits where the centroid and width of the redder Gaussian were fixed to the 
best-fitting values for the full sample. The orange dashed curves are single-Gaussian fits. 



In Figure [2l we also show a single- Gaussian fit to the 
Ca H&K feature. Besides being a poor representation of the 
data, the centroid of the Gaussian is consistently intermedi- 
ate to the two components. Usually, one wants to measure 
the photospheric velocity for a given feature. With that goal, 
the single Gaussian clearly fails. A single Gaussian, by mea- 
suring something intermediate to the two components, mea- 
sures nothing physical. Furthermore, the centroid of the sin- 
gle Gaussian is significantly affected by the blue component. 
The single Gaussian fits for the subsamples indicate that the 
low-stretch SNe have significantly lower velocities than the 
high-stretch SNe. However, the double-Gaussian fits show 
that this is not the case for the photospheric component. 

To investigate the importance of the blue component to 
the measured ucaH&K from these two methods, we created 
artificial, but realistic, line profile s. In Figure O we again 
show the median spectrum from the lMld sample. We created 
a double-Gaussian line profile to mimic the profile of the 
median spectrum. We then varied the height of the bluer 
Gaussian, but left all other parameters fixed. We display 
several example line profiles in Figure O Visually, all of these 
line profiles appear physically possible and represented in 
nature. The full sample of line profiles vary from having no 
blue component to having a blue component that is about 
twice as strong as the red component. 

We fit single Gaussians to all artificial line profiles. We 
display a subset of these fits in Figure [3] (those that match 
the subset of profiles displayed). As expected, the stronger 
the blue component, the bluer the centroid of the Gaussian. 
In FigurelH we show the measured vca. h&k from these Gaus- 
sian fits. Over the range we explore (from no blue compo- 
nent to a blue component that is twice as strong as the red 
component), the measured «caH&K changes by more than 
5000 kms"'^. Even when the blue component is about a 
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Figure 3. Ca H&K li ne pro files. The black curve is the median 
spectrum from the full lM12l sample. The solid lines are artificial 
line profiles created from two Gaussians where only the height of 
the blue component varies. The dotted lines are single-Gaussian 
fits to the artificial profiles. 



third as strong as the red component, the measured «caH&K 
is ~1000 kms~^ different from the true wcaH&K- 

We also measured the wavelength of maximum absorp- 
tion. This wavelength is associated with the blue component 
when it is stronger and quickly transition its association to 
the red component as the blue component becomes weaker. 
The measured uca h&k for our artificial line profiles is shown 
in Figure [l] Although this method fails dramatically for 
strong blue components, the measured ncaH&K is relatively 
constant for line ratios less than one, with all measured ve- 
locities <1000 kms"'^ for all such cases. There is a slight 
bias (~150 kms~^) for these measurements, some of which 
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Figure 4. Measured velocity for artificial line profiles. The hori- 
zontal black lines represent the velocity of the two components (as 
measured from their centroid and assuming a rest wavelength of 
3945 A), with the lower and higher velocity components labelled 
"Photospheric Calcium" and '"High Velocity' Calcium," respec- 
tively. The black crosses represent the measured vcaH&K from 
a single Gaussian to fit the profiles. The blue X's represent the 
measured "UCaH&K from the wavelength of maximum absorption. 



can be explained by increasing flux of the pseudo-continuum 
with wavelength. Correcting for the pseudo-continuum re- 
moves much of the bias, with the remaining bias related to 
the strength of the blue component. 

For cases where the red component is stronger than 
the blue component, measuring the wavelength of maximum 
absorption is significantly better at measuring the photo- 
spheric velocity than using a Gaussian fit to the full profile. 
In this regime, the wavelength of maximum absorption is 
only minimally affected by the strength of the blue com- 
ponent, while the Gaussian fit is significantly affected. In 
the regime of having a stronger blue component, the wave- 
length of maximum absorption fails. However, in this regime, 
the Gaussian fit also fails, producing unphysical and signif- 
icantly biased results. 

Using the lFolev etHI l|201ll ) method of culling uca h&k 
measurements that are not representative of the photo- 
spheric velocity, one should have reliable iiCaH&K measure- 
ments, but will necessarily have an incomplete sample. A po- 
tential way to avoid this bias would be to perform a double- 
Gaussian fit. We have not investigated how this method per- 
forms with noisy data. 



3 SYNOW MODELS 

To further understand the nature of the Ca Hfc K feature, we 
use t he SN spectrum-synthesis code SYNOW (jPisher et al.l 
[1993) to create simple SYNOW spectral models. We specif- 
ically use these models to test how temperature can affect 
the profile and look for trends between the Ca H&K profile 
shape and other spectral features. Although SYNOW has 
a simple, parametric approach to creating synthetic spec- 
tra, it can provide insight on basic trends in SN SEDs. To 
generate a synthetic spectrum, one inputs a blackbody tem- 
perature (Tbb), a photospheric velocity (iiph), and for each 
involved ion, an optical depth at a reference line, an ex- 



citation temperature (Tcxc), the maximum velocity of the 
opacity distribution (iimax), and a velocity scale {ve). This 
last variable assumes that the optical depth declines expo- 
nentially for velocities above «ph with an e-folding scale of 
Ve- The strengths of the lines for each ion are determined by 
oscillator strengths and the approximation of a Boltzmann 
distribution of the lower-level populations with a tempera- 
ture of Texc. 

We produced models consisting of only Call and with 
only Sill and Call to isolate their affect on the profile of 
the Ca H&K feature. For aU models, we set Tbb = 10000 K, 
?;ph = 10000 kms-\ u^ax = 80000 kms'S and = 3000 
(for Sill) and 2000 kms"^ (for Call). We chose r = 5 and 4 
for Sill and Call, respectively. These parameters were cho- 
sen such that when Tcxc = 10000 K, the model Ca H&K line 
profile was visually similar to that of the median spectrum 
of the iMld sample. Keeping all other parameters fixed, we 
varied Texc from 5000 to 20000 K. A subset of the models 
spanning this range are presented in Figure [5] 

As seen in Figure [5] the inclusion of Si II dramatically 
changes the Ca H&K profile shape, making it stronger, 
broader, and bluer. Although the Sill A3858 feature may 
be stronger in the models than in real SN spectra, the Sill 
A6355 and the Ca H&K features appear to have reasonable 
strengths. 

There is a clear spectral progression as the tempera- 
ture changes. We note that for SYNOW, Tbb only changes 
the continuum shape of the models and does not affect the 
strength of features. Since SYNOW uses Ca H&K as the ref- 
erence calcium line, the strength of the Ca H&K absorption 
by definition does not change much with Tcxc, and the entire 
calcium spectrum does not change much over the temper- 
atures probed. Meanwhile the Sill spectrum changes sig- 
nificantly with varying Tcxc. The strength of the Ca H&K 
absorption within the Ca H&K feature (i.e., the strength 
of the red component) does change slightly with Texc be- 
cause of the strength of the Sill A3858 emission changing 
the apparent Ca H&K absorption. 

In the red, there is the expected change in the ra- 
tio of the Sill A5972 and A6355 lines. This ratio, 7^(Si), 
is highly correlated with luminosity and light-curve shape 
(|Nugent et al.l Il995h . As the Sill A5972 feature becomes 
stronger, the Sill A3858 feature becomes weaker. For the 
SYNOW models, 7?.(Si) increases with increasing Texc, while 
for SNe la, 7^(Si) incre ases with decreasing T ; this has been 
previously noted (e.g.. iBongard et al.ll2008l ). and is likely 
the result of not simultaneously changing the opacity with 
Texc and/or non-local thermodynamic equilibrium effects. 
However, the Call spectrum does not change significantly 
with Tcxc and other model spectra show t he same relation 
between Sill A3858 and Sill A5972 (e.g., iKasen & Plewal 
[2007; Blondin ct al. 2012a). We therefore consider the qual- 
itative changes in the spectra to be correct, although the 
corresponding temperatures may not be. All models show 
that the strength of Sill A3858 and Sill A5972 are anti- 
correlated; we will use this relation as the primary model 
prediction. We will later use 7?,(Si) as a proxy for light-curve 
shape. 

The excitation energy for the various Si II lines also ex- 
plain the correlations between the various Sill features. The 
Sill A3858, Sill A5972, and Sill A6355 features have excita- 
tion energies of 6.9, 10.0, and 8.1 eV, respectively. Because 
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Figure 5. SYNOW model spectra. The dashed and sohd curves represent models including only Call and both Sill and Call, respec- 
tively. The models only vary in their excitation temperature, which is labelled. The midd le and left panels show detailed views of the 
Ca H&K feature and the redder Sill complex, respectively. The median spectrum from the lMl j sample is shown in the middle panel to 
demonstrate that the Tcxc = 10,000 K model has a similar Ca H&K profile shape. 



the Sill A3858 and Sill A5972 features have very different 
excitation energies and Si II A6355 has an excitation energy 
intermediate to the other two features, the strengths of the 
Sill A3858 and Sill A5972 features should change in oppo- 
site directions with changing temperature. This also explains 
the SYNOW results since SYNOW fixes the strength of the 
reference feature, Sill A6355. 

In the middle and right-hand panels of Figure O we 
show the Ca H&K feature and redder Si II complex in detail. 
Again, it is clear that both 7^(Si) and the strength of the 
Sill A3858 feature change in the way described above. 

A SN photosphere is, of course, more complicated than 
the simple SYNOW model. Specifically, as the temperature 
changes over the relevant range, the ionization of silicon 
(specifically the amount of singly and doubly ionized silicon) 
changes. Additionally, certain features may be saturated 
(and possibly for only certain temperatures). Specifically, 
it is thought that Sill A6355 is usually saturated. However, 
it seems unlikely that the other Sill features are saturated. 
Therefore, these additional complexities should not change 
our interpretation that the strengths of the Sill A3858 and 
Si II A5972 features are anti-correlated and change with tem- 
perature. 

We fit the Sill A4130, A5972, and A6355 features in 
each model spectrum with single Gaussians. Although the 
line profiles are not exactly Gaussian, the fits are reasonable 
approximations of the data, and the process is similar to 
what is done in practice. We also fit the Ca H&K feature 
with both a single Gaussian and a double Gaussian. We show 
the measured velocity in Figure [6] 

The measured Sill A4130, and A6355 velocities differ 



by at most 530 and 220 kms~^ over the entire temperature 
range, respectively. At the lowest Texc, the Sill A5972 fea- 
ture is not strong enough to measure a reliable velocity, but 
for the other temperatures, it differs by at most 440 kms~^. 
These differences are encouraging since the photospheric ve- 
locities did not change. 

Fitting two Gaussians to the Ca H&K feature, which 
should be better at recovering the true velocity (see Sec- 
tion[2}, we see that the red component, corresponding to Ca 
H&K, has a measured velocity range similarly small to that 
of the Sill features noted above. Specifically, the maximum 
difference of measured Ca H&K velocities over all tempera- 
tures probed is only 290 kms~^. However, the measured 
velocity for the Sill A3858 feature changes significantly 
with temperature. Over the full temperature range probed, 
the measured Sill A3858 velocity ranges from —15,550 to 
— 19,380 kms~^ - a difference of 3820 kms~^, or roughly 
an order of magnitude greater than that of the other fea- 
tures. 

A single- Gaussian fit performs even worse. Because of 
the changing Sill A3858 velocity and its varying strength, 
a single- Gaussian fit to the Ca H&K feature results in 
a fCaH&K range of —13,250 to —17,670 kms~^ over our 
chosen temperature range for a maximum difference of 
4420 kms~^. Since the ucaH&K measured using a single 
Gaussian can have dramatic differences even when there is 
no change in physical velocities, this is even more reason to 
avoid this technique. 

Figure[n]also shows the Sill A5972 to Sill A6355 (7^(Si)) 
and Sin A3858 to Ca H&K ratios, which we will call the 
Si/Ca ratio. The range for 7?.(Si), from effectively zero (when 
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Figure 6. Top Panel: Measured velocity as a function of excita- 
tion temperature for SYNOW model spectra. The red dotted lines 
represent the velocities measured v^ith a single Gaussian for the 
Si II A4130 and A6355 lines. The black crosses represent velocities 
measured with a single-Gaussian fit to the Ca H&K feature. The 
blue X's represent velocities measured with a double-Gaussian fit 
to the Ca H&K feature, where the Si II A3858 velocity is measured 
assuming its rest wavelength is 3945 A. Bottom Panel: The 7?,(Si) 
(red dotted line) and Si/Ca (blue X's) line ratios as a function of 
excitation temperature for the SYNOW model spectra. 



Sill A5972 is difficult to discern) to ~0.3, is approximately 
the range seen by all SNe la except SN 1991bK-lik e ob- 
jects (e.g.. lBlondin et al.|[2012bl : ISilverman et al.ll2012l) . The 
Ca/Si ratio has a range of 1.2 to 2.3. The ratio is affected by 
both the strength of the Si II A3858 absorption and the Si II 
A3858 emission, which fiUs in some of the Ca H&K absorp- 
tion. As noted above, the strength of the Sill A3858 feature 
can have a large affect on the Ca H&K line profile, and even 
dominates for many temperatures. 

Since Sill A5972 and A6355 do not show significant ve- 
locity differences with temperature, are relatively free of 
contamination from other species, and 7?.(Si) is a good in- 
dication of light-curve shape, we can use it as a proxy for 
light-curve shape for our models. Figure [7| shows «ca h&k 
measured with a single Gaussian as a function of 7?.(Si). 
The measured iicaH&K decreases in amplitude with increas- 
ing 7?.(Si), which corresponds to decreasing stretch and lu- 
minosity. 

Converting stretch to 7^ (Si), w e can plot the lM12l mea- 
surements on Figure [T] The IM12| spectra do not cover the 
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Figure 7. Velocities measured with a single-Gaussian fit to the 
Ca H&K fea ture a s a function of 7e(Si) for the SYNOW model 
spectra. The |M12| data are also plotted in blue where the mea- 
sured stretch is converted to 7?.(Si) with the top axis showing the 
scaling. 



redder Sill featu res, so a direct measurement could not be 
made. The IM12| values, which use a single Gaussian to fit 
the Ca H&K feature, span a similar range of wcaH&K and 
inferred 7?,(Si) as the models, with the model trend going 
through the middle of the data values. The claime d tren d of 
«CaH&K with light-curve shape is clear using the IMI j val- 
ues. However, this trend is similar to the trend generated 
by simply changing the temperature in the SYNOW models. 
We emphasize that the true wcaH&K is fixed for all models 
and the wcaH&K measured using a double- Gaussian fit varies 
only slightly over all models. The trend shown in Figure [7] 
is solely the result of the method for measuring the velocity. 
The underlying physical effect is the changing strength of 
Sill A3858 with temperature, and thus light-curve shape. 

There will undoubtedly be a true range of velocities in 
the data. We have not explored how the SYNOW models 
change when varying se veral parameters, but it is clearly 
possible to reproduce the lM12| trend through a combination 
of single Gaussian fitting, an inherent velocity range, and re- 
lations between the Ca H&K profile shape with temperature 
and velocity. 



4 THE M12 SAMPLE 

With the insigh ts of the above analysis, we re-analyse the 
individual 'M1^ spectra. In Section O we showed that fitting 
the Ca H&K profile with a single Gaussian results in an 
imprecise, biased, and unphysical measurement of the pho- 
tosph eric velocity. Instead, we fit the Ca H&K profiles of the 
IM12| spectra with two Gaussians. As suggested above, this 
method should provide relatively unbiased measurements of 
the Si II A3858 feature and the Ca H&K feature. 

We provide the best-fitting velocities for Si II A3858 as- 
suming that it is "HV" Ca H&K and Ca H&K in Table [T] 
We also provide the Si/Ca ratio in Table [T] 

Figure |8] compares the velocities measu red with the 
double-Gaussian fit to those reported bv lM12l . The vcaH&K 
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Figure 8. IM12| Ca H&K velocities measured using a single- 
Gaussian fit to the Ca H&K feature compared to velocities mea- 
sured with a double-Gaussian fit to the Ca H&K feature for the 
IM12| sample. The black and blue points represent the Ca H&K 
and Si II A3858 (assuming a rest wavelength of 3945 A) velocities 
from the double-Gaussian fit. The solid lines represent the best- 
fitting linear relationships for the data, corresponding to 1.5 and 
3.1-(T results for Ca H&K and Sill A3858, respectively. 
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Figure 9. Velocities measured with a single or double-Gau ssian 
fit to the Ca H&K feature as a function of stretch for the 
sample. The single- Gaussian measurements are taken from 



M12 



M12 



and represented by the red data. The double-Gaussian measure- 
ments are shown as black and blue points for the Ca H&K and 
Sill A3858 (assuming a rest wavelength of 3945 A) components, 
respectively. The solid lines represent the best-fitting linear rela- 
tionshi ps for the data, corresponding to 2.3, 2.3, and lA-a results 
for the lMlSl Si II A3858, and Ca H&K velocities, respectively. 



is systematically lower tlian the lM12l value for the same SNe. 
Similarly, tlie Sill A3858 (whe n trea ted as HV Ca H&K) is 
systematically higher than the lMld value for the same SNe. 
Perfo rming a Bayes ian Monte-Carlo lin ear regressio n on the 
lM12l -Sin and lM12l -Ca H&K data sets (|Kellvl(2007h . we find 
that 99.8 and 86.8 per cent of the realizations have positive 
slopes, respectively, corresponding to 3.1-a and 1.5-cr results, 
respectively. That is, ther e is s ignificant evidence for a fin- 
ear relation between the IM12| vcaH&K measurements and 
the Si II A38 58 ve locity, but no evidence for a linear relation 
between the lMld ucaH&K measurements and the velocity of 
the red component corresponding to absorption from pho- 
tospheric calcium. 

As expected from the results of Section |3l the IM12| 
fCaH&K measurements are intermediate to the Sill A3858 
and Ca H&K velocities, more closely track the Sill A3858 
velocity than the Ca H&K velocity, and are systematically 
biased measurements of wcaH&K- 

For much of the analysis performed bv lMlj . they cor- 
rected their measured vcslHScK to a maximum-light value, 
iiCaH&K) using a single velocity gradient for all objects de- 
rived from the wcaH&K and eflective phase measurements 
for their sample. Using a large sample where many ob- 
jects had multiple spectra obtained near maximum light, 
iFolev et all (120111 ) showed that 

"^CaH&K and the velocity 
gradient were highly correlated with higher-velocity SNe 
also having higher- velocity gradients. Taking into account 
this relation, they produced an equation w hich e stimates 
Wca H&K given uca H&K aud phase. Using the IM12| velocity 
gradient results in differences between vcaH&K and VcaH&K 
that can be as large as 1150 kms"'^, and the median dif- 
ferenc e is 460 kms~^. However, if one uses the lFolev et al.l 
relation to correct the velocity measurements to have 
a common phase of 2.7 d (the median of the sample), then 



the deviation between that value and ucaH&K is at most 
560 km s~^, with a median absolute deviation of 120 kms~^, 
both of which are smal ler than the typical uncertainty of 
670 kms"^ for the lM12l 

'^CaH&K measurements. The combi- 
nation of using a single velocity gradient and extrapolating 
to maximum light (only one spectrum in the sample has a 
phase before maximum brightness) introduces unnecessary 
additional uncertainty. Instead in all further analysis, we 
use the raw «CaH&K measurements, but add an additional 
120 kms~^ uncertainty in quadrature to the reported un- 
certainty. 

Using our measurements of ucaH&K, we can re-examine 
the IM12| claim that wcaH&K is correlated with light-curve 
shap e. In Figure |9l we show the Sill A3858, Ca H&K, and 
IM12| velocity measurements as a function of stretch. This 
figure is similar to Figure [T] 

Pe rform ing a Bayesian Monte-Carlo linear regression 
on the iMll the Sill A3858, and Ca H&K velocity mea- 
surements, we find that 97.7, 97.8, and 87.6 per cent of the 
realizations have positive slopes, respectively, corresponding 
to 2.3-a, 2.3-CT, and 1.5-cr res ults, r espectively. We therefore 
find mild evidence that the IM12| and Si II A3858 velocity 
measurements are linearly related to stretch. We find no sta- 
tisti cal ev idence that wcaH&K is linearly related to stretch. 

IM12| found a 3.4-cr linear relation between their ^caH&K 
measurements and stretch. Above, our measured significance 
is much lower. This is partly because in the above calcula- 
tion, we did not in clude SNe 2011by and 2011fe since we do 
not have the lM12| spectra. Including the reported values for 
these SNe, there is only a minor change in the significance, 
changing the percentage of realizations with positive slopes 
to be 98.7 per cent, which is a 2.5-cr result. The other differ- 
ence is that above we examine wcaH&K instead of WcaH&K- 
Performing the same analysis as IM12| (using 
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eluding SNe 2011by and 2011fe), we find that 99.4 per eent of 
the realizations have positive slopes, which is a 2.8-(j result. 
The difference in significance is likely i n the subtleties of fit- 
ting a line. Thi s pract i ce is n ot trivial (|Hogg. Bow. &: Land 
|2010| ). but the iKellvl (|2007l ') method is generally a better 
choice than most options. 

We also performed a Kolmogorov-Smirnov (KS) test, 
splitting the samples by a stretch of 1.01. This is not an 
ideal test since there are several SNe with stretches consis- 
tent with 1.01 (and therefore could be in either group) and 
uncertainty in the velocity can also change the overall distri- 
bution, but it can provide an indication of a difference. The 
KS test r esuhed in p values of 0.0014, 0.0036, and 0.080 for 
the |M12| . the Sill A3858, and Ca H&K velocity measure- 
ments, respectively. These tests indicate that the low and 
high-stretch subsamples ha ve diff erent parent populations 
for both the Sill A3858 and lM12l velocities. However, there 
is no statistical evidence that the low/high-stretch subsam- 
ples have different parent ucaH&K distributions. 

Although there is only margi nal ev idence that there is 
a linear relationship between the IM12| measurements and 
stretch, we find a similar significance of a relationship be- 
tween Sill A3858 velocity and stretch. Since the vcaH&K 
shows no evidence for a correlation with stretch and the lM12l 
measurements are correlated with the Sill A3858 velocity 
(and at most weakly correlated with the Ca H&K velocity), 
the p hysical relationship underlying the result identified by 
IM12| is likely the correlation between Sill A3858 velocity 
and stretch. From the SYNOW models, this relation is un- 
derstood as a temperature effect (and not a real difference 
in photospheric velocity). 

Next, we examine the Si/Ca ratio. We show the Si/Ca 
ratio as a function of stretch for the lM12l sample in Figure fTUl 
Performing a Bayesian Monte-Carlo linear regression on the 
data, we find that 82.4 per cent of the realizations have pos- 
itive slopes, corresponding to a lA-a result. Although there 
is no evidence for a li near relationship between stretch and 
the Si/Ca ratio in the lM12l data, there is a slight correlation 
with higher stretch SNe having larger Si/Ca ratios. This is 
the same general trend expected from the SYNOW models, 
and future investigations should determine if such a trend 
exists. 

Finally, we compare the Si/Ca ratio to our measured 
velocities (Figure [TTjl . There are no obvious trends (1.0 and 
1.7-0-) between the Si/Ca ratio and the Sill A3858 or Ca 
H&K velocities. However , ther e is a moderate trend between 
the Si/C a ratio and the IMI j measurements (2.7-cr), where 
the lM12l velocities increase with increasing Si/Ca ratio. One 
sh ould e xpect that the single-Gaussian method (as employed 
bv lM12t ) should be intermediate to the Sill A3858 and Ca 
H&K velocities. The velocity should be closer to the Ca 
H&K velocity when the Sill A3858 feature is weak (small 
Si/Ca ratio) and closer to the Sill A3858 velocity when the 
feature is strong (large Si/Ca ratio). 

The behavior seen in the data is reproduced to some 
extent by the SYNOW models. Figure [11] also shows the 
SYNOW model Si/Ca ratio compared to the Sill A3858, Ca 
H&K, and single-Gaussian velocities. The Ca H&K veloc- 
ity is relatively flat for the SYNOW models, while the Sill 
A3858 and single-Gaussian velocities increase with increas- 
ing Si/Ca ratio. The Ca H&K and single-Gaussian trends 
are similar in the data and models, but the data have higher 
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Figure 10. Si/Ca ratio as a function of stretch for the lMl j sam- 
plc. The solid line represents the best-fitting linear relationship 
for the data, corresponding to a 1.4-cr result. 
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Figure 11. Velocities measured with a single or double-Gaussian 
fit to the Ca H&K feature as a function of the Si/Ca ratio for the 
M12 sample. The single- Gaussian measurements are taken from 
MI2I and represented by the red data. The double- Gaussian mea- 
surements are shown as black and blue points for the Ca H&K and 
Sill A3858 (assuming a rest wavelength of 3945 A) components, 
respectively. The solid lines represent the best-fitting linear rela- 
tionships for the data, corresponding to 2.7, 1.0, and 1.7-(t results 
for the M12, Sill A3858, and Ca H&K velocities, respectively. The 
dotted lines represent the best-fitting linear relationships for the 
data with a Si/Ca ratio of >1. The red, black, and blue circles 
represent the single-Gaussian and double-Gaussian fit velocities 
for the SYNOW model spectra. 



velocities than the data (by about 700 and 2200 kms~^, re- 
spectively). There is no obvious trend in the Sill A3858 data, 
contrary to what is seen in the models. 

The models are restricted to Si/Ca ratios of >1. If we fit 
the data with the same restriction, the fits are more similar 
to the slopes seen in the models. Specifically, the trend with 
Ca H&K is flatter and the trend with Si H A3858 is stronger. 
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although neither trend is statistically significant. Perhaps a 
simple linear relation is not sufficient to describe the trend 
between the Si/Ca ratio and Sill A3858 velocity. 



5 THE CFA SAMPLE 



Although the lM12l sample was used to identify some trends 
which were investigated above, it is limited in size. To 
test additional trend s, we use the CfA spectral sample 
iBlondin et al.ll2012bl ). Over the last two decades, the CfA 
SN Program has observed hundreds of SNe la, m ostly with 
the FAST spectrograph (|Fabricant et al.l 1 19981 ) mounted 
on the 1.5 m telescope at the F. L. Whipple Observa- 
tory. 'rhe_d^ta^Jiwe_beenj;^^ man- 
ner l|Matheson et al.ll2008l : [siondin et al.ll2012bl '). producing 
well-calibrated spectra. These spectra often cover Ca H&K 
and always cover Si II A6355, which the lM12l spectra do not 
cover. 

For SNe la in the sample with a measured time of 
maximum brightness f rom light curves, vsi i i and vqi, h&k 
have been measured (|Blondin et al.l l2012lj ). Briefly, this 
is achieved by first generating a smoot hed spectrum 
using an inverse-variance Gaussian filter ijBlondin et al.l 
l2006h . and the wavelength of maximum absorption in the 
smoothed spectrum is used to determine the velocity (see 
iBlondin et al.l2012bl for details). T he measurem e nts fo r each 
spectrum have been reported by iFolev et all (I2OIII). and 
measurements in all cases were obtained bv iBliandin et al] 
(|2012bl ) 

The CfA sample contains 1630 vsi 11 and 1192 vca h&k 
measurements for 255 and 192 SNe la, respectively. 

The velocity of each absorption minimum in the Call 
H&K fea ture of the smoothe d spectra is automatically 
recorded (|Blondin et al.|[2012bl ). We only examine spectra 
with one or two minima. If two minima are found, the 
higher/lower velocities are classified as "blue" /"red." If only 
one minimum is found, it is categorized as the red or lower- 
velo city component. 

iFolev et all l|201ll ) noted that comparing all nca h&k 
measurements to their corresponding vsi 11 measurements, 
there were two distinct "clouds" corresponding to a lower 
and higher velocity relative to vsi 11. The higher- velocity 
cloud typically corresponds to the blue velocity component, 
although there are some red measurements in that cloud. 
The red measurements in the blue cloud typically have indi- 
cations of a lower-velocity component, such as a red shoulder 
in the line profile, and it was assumed that they likely corre- 
sponded to measurements which were physically similar to 
the blue measurements and were simply misclassified. 

In Figure 1121 we show the subset of CfA measure ments 
of spectra with —1 ^ t ^ 4.5 d (chosen to match the iMlsi 
sample). This subset also s hows the distinct blue/red clouds. 
We u sed the method of IWilliams. Bureau. &: Cappellaril 
l|20ld l to fit a single slope, but separate offsets to the two 
clouds. As a result of that fitting, there is a natural dividing 
line between the two clouds, and we used this line to produce 
cleaner subsamples. We removed every blue measurement in 
the red cloud, since they may be errant measurements. We 
also reassigned every red measurement in the blue cloud as 
a "blue" measurement because of the reasons listed above. 



This full process is shown graphically in the three panels of 
Figure [H 

With these clean subsets, we have a reasonable estimate 
of the velocities for the blue and red components of the Ca 
H&K feature. Since some SNe in the CfA sample have multi- 
ple spectra in the chosen phase range, we created samples for 
each velocity group where there is one measurement per SN. 
For each SN, we chose the mea surem ent closest to a phase 
of t — 2.7 d, the median of the IMI j sample. This resulted 
in samples of 66 and 67 SNe la with blue and red mea- 
surements (approximately one-t hird o f the full CfA sample 
and about 5 times as large as the lM12l sample), respectively. 
We present those measurements as a function of light-curve 
shape (specifically, Ami5(_B)) in Figure [131 There is no sig- 
nificant linear relation between light-curve shape and the 
individual velocity components. In fact, the stronger rela- 
tionship of the red velocity component, which is the best 
representation of the Ca H&K photospheric velocity, is a 
1.3a result in the opposite direction than the lMl j relation 
(i.e., higher velocity for slower-declining SNe la). Because of 
this opposin g tren d, the CfA data are significantly inconsis- 
tent with the M12l relation. Ho weve r, this result is consistent 
with that of lFolev et all (|201ll ') and lFolevI (|2012l ') for both Ca 
H&K and Sill A6355. 

In addition to the arguments detailed above, there is 
additional evidence in the CfA data that suggest that the 
blue component of the Ca H&K feature is the result of Si II 
A3858. If the blue component is predominantly from HV Ca 
H&K, then one would not expect a particularly high correla- 
tion between its velocity and that of Sill A6355. That is, the 
velocity of a HV calcium component could be independent 
of the photospheric silicon velocity. However, there is a rea- 
sonable correlation (correlation coefficient of 0.54) between 
the two. On the other hand, the velocities of the red and blue 
components are barely correlated (correlation coefficient of 
0.27). Therefore, the velocity of the blue component has a 
larger association with the photospheric velocity of silicon 
than the photospheric velocity of calcium. 

Perhaps the best evidence that the blue component is 
from Sill A3858 absorption is presented in Figure [12] In the 
right panel, we plot the relation between vsm and ucaH&K 
in the scenario where the wcaH&K measurement is from a 
misidentified Si II A3858 feature at the same velocity as Si II 
A6355. The line goes directly through the blue cloud, indi- 
cating that the blue component of the Ca H&K feature has a 
velocity consistent with that of Si II A6355 if it is formed by 
Si II A3858 absorption. In other words, the blue component 
is at the waveleng th one expects by blue shifting 3858 A by 
usiii. Addtionallv. iBlondin et aU (|2012bl ) presented several 
examples of SNe where the blue component was consistent 
with Sill A3858 at the Sill A6355 velocity, while the red 
component was consistent with the velocity of the Ca NIR 
triplet. 

From the large CfA sample, we showed additional ev- 
idence that iiQaH&K does not correlate with light-curve 
shape. The velocity of the blue component is correlated with 
the photospheric silicon velocity (as measured by Si II A6355) 
and relatively uncorrelated with the photospheric calcium 
velocity. In addition to being correlated with photospheric 
silicon velocity, the CfA data show that the velocity of the 
blue component matches the expected photospheric velocity 
of Sill A3858. 
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Figure 12. Velocities from the Ca H&K profile (both blue and red components), assuming a rest wavelength of 3945 A, as a function 
of Sill A6355 velocity for the CfA sample with a phase range of —1 t ^5 4.5 d. The middle panel plots the data as blue and red for the 
blue and red compon ents, respectively. The dashed line represents the line of separation for the two clouds from the fitting method of 
IWilliams et al.l hOl(t) . The right panel removes all measurements from blue components in the red cloud (as determined by the dividing 
line) and assigns all "red" velocities in the blue cloud as "blue" points under the assumption that they were originally misclassified. The 
dotted line represents spectra that would have the same Ca H&K and Sill A6355 velocities. The solid line represents spectra that would 
have the same Ca H&K and Si II A6355 velocities if Ca H&K were actually from Si II A3858 absorption at the same velocity as Si II 
A6355. 



6 ADDITIONAL MODELING 

From the abo ve an alysis of the SYNOW models, the com- 
parison of the IM12| sample to the SYNOW models, and an 
examination of the CfA sample, there is significant evidence 
that the blue component of the Ca H&K feature is predom- 
inantly from Sill A3858. However, additional confidence in 
this claim can be obtained by modeling specific SNe. 

In this section, we examine the two possible scenarios 
for the blue component of the Ca H&K feature (either HV 
calcium or Sill A3858) for two test cases: SN 2011fe and 
SN 2010ae. 



6.1 SN 2011fe 

SN 2011fe, which occurred in M 101 and was the bright- 
est SN la in 40 years, has been incredibly well ob- 



served and extensively studied (e.g., 


Nueent et al. 


20 111; 


Brown et al.ll2012l; IChomiuk et al.il2012 


: iHoresh et al.l 


2OI2I: 


MarEutti et al.l 20121: 


Matheson et al.l 


2OI2I: Parrent et all 


20121: IShaDoee et al.l 


20131). Here we examine a single 



maximum-light spectrum of SN 2011fe, obtained by HST 
using the STIS spectrograph (Program GO-12298; PI El- 
lis). The spectra were obtained on 2011 September 10 be- 
tween 09:51 and 11:14 UT, correspond ing to t = 0.0 d rel- 
ative to B-band maximum brightness fMlz). The observa- 
tions were obtained with three different gratings, all with 
the 52" X 0."2 slit. Two exposures were obtained for each 
of the CCD/G230LB, CCD/G430L, and CCD/G750L se- 
tups with individual exposure times of 530, 80, and 80 s, 
respectively. The three setups yield a combined wavelength 
range of 1665 - 10,245 A. The data were reduced using the 



standard HST Space Telescope Science Data Analysis Sys- 
tem (STSDAS) routines to bias subtract, fiat-field, extract, 
wavelength-calibrate, and fiux-calibrate each SN spectrum. 

We present the spectrum in Figure 1141 We note that 
IM12| presented a spectrum for SN 2011fe from a different 
phase and which only covered ~2900 - 5700 A. This is the 
first publication of these data. This is also only the second 
published maximum-light SN la spectrum to probe below 
~2500 A (the first being of SN 201 liv; Foley et al. 2012c). 

The SN 2011fe spectrum is of extreme high quality, in- 
cluding in the UV. Because of its quality and wavelength 
coverage, we can produce a reasonable SYNOW model. We 
have made two attempts at fitting the SN 2011fe spectrum 
using SYNOW models. The first assumes that the blue com- 
ponent of the Ca H&K feature is from Si II A3858; the second 
assumes that it is caused by HV Ca H&K. We present our 
models in Figure [14] and model parameters in Table [S] 

When generating these models, we first attempted to fit 
the full spectrum with a limited number of species. These 
models are not optimized to fit the entire spectrum; because 
of potential effects other species could have on the spectral 
features of interest, we wanted a first-order model of the full 
spectrum. We then either added HV Call or adjusted the 
Sill temperature to match the blue component of the Ca 
H&K feature. We allow the opacity and density structures 
for Call, Sill, HV Call, and Nal to vary, but all other 
species remain the same. 

For the HV calcium model, we adjusted the Sill tem- 
perature to an extreme value that still fits the Sill A5972 
feature. In this model, we do not include any Nal, and there- 
fore Na D does not contribute at all to this feature. As a 
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Figure 14. HST spectrum of SN 2011fe at t = 0.0 d (black curve). The blue and red curves are "best-fitting" SYNOW model spectra 
where there was an attempt to simultaneously match the Ca H&K profile and Sill A5972 and A6355 lines using Sill/HV Call to match 
the blue component of the Ca H&K feature, respectively, and differences in the strength of Na D to match the Sill A5972 feature. 
The dotted red line represents the model represented by the solid red line except without a HV Call component. The dotted blue line 
represents the model represented by the solid blue line, except with the Call opacity and density structure matched to that of the model 
represented by the solid red line. 



result, the Sill A3858 is about as weak as possible, and the 
HV Ca H&K is essentially as strong as possible. 

For the Sill A3858 model, we adjust the Sill tempera- 
ture to an extreme value to match the blue feature in the 
Ca H&K feature. We then add Nal to match the strength 
of the feature near 5800 A. 

These model s differ in some ways from those presented 
by iParrent et al.l l|2012l ) for their optical-only maximum- 
light SN 2011fe spectrum. Most differences are related to 
matching the UV region, which requires adding Coll and 
CrII. Interestingly, adding these features reduces the need 
to include Fell in the SYNOW model (although we cannot 
definitively say that it is not in the spectrum). Addition- 
ally, we are able to b etter model the Ca H&K feature than 
iParrent et all (|2012l ) because of the additional data blue- 
ward of the feature. 

Examining the SYNOW models in detail, particularly 



near the Ca H&K feature, the redder Sill features, and the 
Ca NIR triplet (see lower panels of Figure I14p , we see that 
the models are very similar. In other words, SYNOW mod- 
eling of SN 2011fe cannot distinguish between our two sce- 
narios; it simply has too many parameters for the data. 



We did not adjust the models to fit the Ca NIR triplet, 
with the hope that we might see signatures of HV Ca. There 
is a feature in the SYNOW model that is coincident with a 
shoulder in the SN 2011fe spectrum. However, we see a sim- 
ilar feature in the Si II A3858 model that is simply the result 
of a slightly different density profile for Call. A full spec- 
tral sequence and /or NIR spectra, which would supply addi- 
tional Sill features, may provide a clear way to distinguish 
the models. 
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Figure 13. Velocity of the blue (top) and red (bottom) com- 
ponents of the Ca H&K profile (assuming a rest wavelength of 
3945 A) vs. Ami^{B) for the CfA sample and a phase range of 
— 1 t 4.5 d. The choice of blue and red velocities were made 
as described in the text and visually represented in Figure 1121 
Each point represents a different SN (66 and 67, respectively). 
There is no significant correlation for either (linear relations are 
0.4(T and l.Scr significant, respectively). Since the red component, 
corresponding to photospheric Ca H&K, has a slight (but insignif- 
icant) trend of higher velocity with faster-declining lig ht cu rves, 
the data are significantly inconsistent with the claim of lMlj that 
Ca H&K velocity decreases with faster-declining light curves. 



6.2 SN 2010ae 

With an inconclusive result from modeling SN 2011fe, we 
now turn to mode ling SN 2010ae. SN 20 10ae is a SN lax 
jFolev et al.ll2 012a') similar to SN 2008ha ()Foley et al.ll2009l . 
I2OI0I : IValenti ct al. 200^. Its spectrum is similar to that of 
a SN la, but with an extremely low ejecta velocity. This in- 
dicates that the ejecta composition, density structure, tem- 
perature, and other aspects of the explosion important for 
producing a particular SED are similar for SN 2010ae and 
SNe la. However, because of the low ejecta velocity, line 
blending is minimal. 

We present a near maxi mum-light spectrum of 
SN 2010ae originally presented by iFolev et~al1 l|2012ah and 
presumed to be obtained near maximum light in Figure [TS] 
This spectrum only covers optical wavelengths. We dered- 
dened the spectrum by E{B — V) = 0.6 mag to roughly 
match the continuum of SN 2011fe and smoothed the spec- 
trum with a inverse-variance weighted Gaussian filter and 
velocity scale of 150 kms~^. 

Perhaps the most important aspect of the SN 2010ae 



spectrum is that the Ca H&K feature is separated into two 
distinct features. 

We then attempted to produce SYNOW model spectra 
in a way similar to what was performed for SN 2011fe. As a 
starting point, we used the SN 2011fe models. We decreased 
'i'phot from 9000 kms~^ to 3000 kms~^. We also reduced 
minimum and maximum velocities for each species. The de- 
tails of the models are presented in Tabled 

We did not change the majority of parameters for the 
model. As a result, the fits are not ideal. In particular, the 
lack of C II results in missing obvious features. Additional 
adjustments would certainly improve the overall fit, but this 
is not necessary for our purpose. However, keeping the model 
similar to that of a SN la (with mostly just adjustments 
to the velocity) reinforces the spectral (and compositional) 
similarities between SNe lax and SNe la. 

Additionally, we changed the opacity of Sill and Call, 
and we changed the density structure of Call. We adjusted 
the opacity of Sill to roughly match the Sill A6355 fea- 
ture. The Call opacity was changed to roughly match the 
NIR triplet. The velocity of the HV calcium and the density 
structure of both the HV and photospheric calcium were 
adjusted to match the Ca H&K feature. 

Ca H&K are offset by 34.8 A, which corresponds to 
2640 kms^^. The velocity difference between the two com- 
ponents will be present even if Ca H&K are blueshifted. 
For most SNe, the ejecta velocities are high enough where 
the two components blend together completely. But for 
SN 2010ae, which has an ejecta velocity similar to this sepa- 
ration, any Ca H&K feature will be roughly twice the width 
of a feature from a single line. For SN 2010ae, the blue com- 
ponent of the Ca H&K feature has a FWHM of 2960 km s~^ 
Therefore, the Ca H&K components can barely fit within 
the width of the feature (with a velocity of ~11200 kms~^, 
about 4 times that of the photospheric velocity), but then 
the line can only be minimally broadened. That is unphysi- 
cal, but if it were the case, then one would expect two com- 
ponents within the blue component, which is not seen. 

The only other choice is to choose a velocity which 
results in either Ca H or Ca K to have a minimum near 
3800 A. Doing this for Ca H results in a velocity of ^ 
13000 kms~^ and a significant absorption feature at ^ 
3760 A, where no such feature exists. When assigning a ve- 
locity of ~ 10000 kms^^ for HV calcium (such that Ca K 
is at ~38QQ A), there is no gap between the blue and red 
components. Neither option reproduces the observed profile 
for SN 2010ae. 

Alternatively, the Sill A3858 model roughly matches 
the spectrum of SN 2010ae. In particular, it reproduces the 
(now unblended) Ca H&K feature. The HV calcium model, 
on the other hand, does not reproduce a key aspect of the Ca 
H&K feature - its unblended nature. It is reasonably certain 
that Si II A3858 causes the absorption of the blue component 
of the normally blended Ca H&K feature for SN 2QlQae. 

Furthermore, removing the HV Ca from the HV Ca 
model does not have two distinct features. It appears nec- 
essary to have a reasonably strong Sill A3858 feature to 
produce the emission between the two components. 

Since iFolev et al.l (|2012al ) showed that SNe lax have 
very similar spectra to SNe la, except with different ve- 
locities, and since the SN 2011fe SYNOW model roughly 
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Figure 15. Same as Figure [141 except for SN 2010ae and its "best-fitting" SYNOW models. The SN 2010ae spectrum has been smoothed 
shghtly and dereddened as described in the text. The sohd and dashed red hnes are for HV calcium where Ca K and Ca H, respectively, 
are matched to the absorption near 3800 A. 



matches the SED of SN 2010ae (with only differences in the 
velocity), one can extrapolate this result to SNe la. 



7 DISCUSSION & CONCLUSIONS 

We have shown through a re-examination of the lMl^ sample, 
a re-examination of the CfA sample, basic SYNOW model- 
ing, and more thorough SYNOW modeling of SNe 2010ae 
and 2011fe that the blue component of the Ca H&K spec- 
tral feature in near-maximum Ught SN la spectra is typically 
fr om Si II A385 8 abso rption. This was also the interpretation 
of I Wang et all (|2003l ). which has spectropolarimetric obser- 
vations of Ca H&K, Sill A6355, and the Ca NIR triplet, 
providing additional weight to this conclusion. Some pre- 
vious claims that the component is the result of HV Ca 
H&K absorption may require re-examination. The Ca NIR 
triplet has shown HV features for some SNe, although it 
is also possible to reproduce some of these features with a 
different (but still smooth) density profile for calcium (see 



Section [6. If) . Therefore, it is still unclear if HV calcium con- 
tributes to the Ca H&K component, how frequently it does, 
and if that contribution is typically blended with Sill A3858. 
The realization that the blue absorption in the Ca H&K pro- 
file is from Sill A3858 for most SNe la has far-reaching im- 
plications for our understanding of SN la progenitor systems 
and explosion models, which have interpreted the prevalence 
of HV calcium as an indication of specific explosion mech- 
anisms and potentially a tracer of the environment of the 
progenitor system. 

Because the Ca H&K profile is a combination of Ca 
H&K and Sill A3858, «caH&K should not be measured by 
fitting the entire Ca H&K feature with a single (Gaussian) 
component. Regardless of the source of the two components, 
we also show that if one does fit the profile with a single 
Gaussian component that the resulting measurements will 
be unphysical, inaccurate, and highly biased. However, be- 
cause of the true nature of the blue component, a single 
Gaussian fit is particularly biased. 

We confirmed the IM12| result that SNe in their sam- 
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pie have different Ca H&K line profiles based on light-curve 
shape. However, the difference is mostly constrained to the 
blue component, with no evidence for a difference in velocity 
or width for the red component. 

We re-examined t he cla im that DcaH&K is corre lated 
with light-curve shape (|M12| ). Using the reported [mT^ mea- 
surements, we do not find a statistically significant linear re- 
lation, but the KS test does indicate different parent popula- 
tions for low/high-stretch subsamples. When using vcaH&K 
measureme nts fr om the red component of the Ca H&K pro- 
file for the IM12| spectra, there is no statistically significant 
trend between vca h&k and light-curve shape. An analysis of 
the CfA sample also showed that there is no correlation be- 
tween ejecta veloc ity and light - curve shap e, confirmin g the 
previous resuhs of iFolev et all (|201ll ) and lFolevI (I2OI2I ). In- 
stead, the unde rlying physical effect driving the relation be- 
tween the lMld measurements and light-curve shape is likely 
the relation between Sill A3858 a nd te mperature. 

This result implies that the IM12| claim that ^caH&K 
does not correlate with host-gala xy m ass is not supported 
by data. Other claims made by IM12| related to WcaH&K) 
including correlations between Vq^ h^^k and the wavelengths 
or velocities of certain features, should also be re-examined. 

From modeling, there is some indication that the Si/Ca 
ratio should be a strong tracer of temperature and an in- 
dicator of light-curve shape, but this is not verified with 
data. There may also be a relatively low correlation be- 
tween ncaH&K and the pseudo-equi valent w i dth o f the Ca 
H&K feature. This may be why Fol ev et all (|201ll ) did not 
find a relation between the pseudo-equivalent width of the 
Ca HfcK feature and t he intr insic c o lour o f SNe la. 

iFolev et all (|201lh and iFolevI (|2012l ) suggested that 
''CaH&K could be useful for measuring the intrinsic colour of 
SNe la. However, this current analysis shows that this ap- 
proach may be limited by the contamination of Sill A3858. 
At the very least, SNe with very high ejecta velocities will 
have a Ca H&K profile that is a blend of Sill A3858 and 
Ca H&K with no distinct components. At that point, one 
should be ci rcumspect o f the d erived velocity. The culling 
technique of iFolev et al.l (|201lh should reduce the number 
of spectra with velocity measurements contaminated by Si II 
A3858, but relatively low signal-to-noise ratio (S/N) spectra, 
galaxy contamination, and other nuisances, may reduce the 
viability of this option. 

There is a proposal to have a low-resoluti on {R f» 75) 
spectrograph on WFIRST (|Green et al.l I2OI2I '). The main 
purpose of the spectrograph for SN science would be spec- 
troscopic classification and redshift determ ination. Simi- 
larly, the SED Machine (|Ben-Ami et al.ll2012l '). is a proposed 
R « 100 spectrograph to classify thousands of low-redshift 
SNe. Another use of these spectrographs could be to measure 
ejecta velocities. Assuming perfect knowledge of the SN red- 
shift, the precision of the ejecta velocity measurement can 
be limited by spectroscopic resolution. 

To test our ability to determine ejecta velocities with 
different resolutions, we show artificial Ca H&K line pro- 
files that contain two components in Figure [161 One cannot 
distinguish the two components of the profile ai R — 50; 
there are ^4 resolution elements in the feature, which is in- 
sufficient for a full six-parameter fit of a double-Gaussian 
fit. Additionally, the two components are separated by ~ 
6000 kms~^, corresponding to i? ~ c/6000 km s~^ ~ 50. 



At R — 75, one can start to see the effect of the two com- 
ponents in some spectra (i.e., ffat bottoms), but the compo- 
nents are still not clearly separate. A resolution of 100 may 
be the minimal amount to clearly see the effects of multi- 
ple components. But considering additional effects such as 
potential [O II] A3727 emission from the host galaxy contam- 
inating the line profile, one might want a higher resolution, 
such as R — 200, where narrow lines should not significantly 
affect the overall profile shape. 

However, we note that Sill A6355 does not suffer these 
same problems, and R = 75 should provide accurate (and 
reasonably precise) measurements of the ejecta velocity. For 
optical spectrographs, one can easily measure vsm to z — 
0.3. With red-sensitive CCDs and good sky subtraction one 
can use optical spectrographs to measure vsm to z ^ 0.6. 
With NIR spectrographs, one can easily measure vsm to 
« 2 (neglecting the faintness of the SNe). 

For the SED Machine, which aims to classify low- 
redshift SNe, it should also be able to measure vsm- The 
proposed spectrograph on WFIRST would have a wave- 
length range of 0.6 - 2 /im, which should cover vsm to 
z « 3, well beyond t he expected redshift range of WFIRST. 
iRodnev et aP (j2012l ) presented an HST observer-frame NIR 
spectrum of a 2 = 1.55 SN la, SN Primo. The spectrum 
has a low S/N and is low-resolution {R « 130). But using 
the method of iBlondin et all |2006l l. we measure vsi 11 — 
— 11200 ± 900 kms~^ at a phase of 6 ± 3 d, corresponding 
to Vsi u = -11700 ± 100 kms"^ This corresponds, us- 
ing the iFolev et all ()201ll ) relations, to {B 

max Vmax ) — 

0.00 ± 0.07 mag. The uncertainty in the velocity measure- 
ment is dominated by the low S/N of the spectrum, but the 
uncertainty in the intrinsic colour is still dominated by the 
uncertainty and scatter in the velocity-colour relation. None 
the less, SN Primo appears to be have a moderate intrinsic 
colour. This shows the potential of using velocity measure- 
ments for SN la cosmology even if the complexities of the 
Ca H&K profile prevents accurate measurements. 

The additional knowledge of the Ca H&K profile pro- 
vided here is a step toward further understanding of the full 
SED of SNe la. SNe lax, which have compositions similar 
to that of SNe la, can be exceedingly useful for determining 
which specific atomic transitions contribute to SN la spec- 
tra. Because of their low ejecta velocities, SNe lax may pro- 
vide additional insight into the specific contributions from 
various lines for blended SN la features. Similarly, additional 
spectropolarimetric observations of SNe la, and particularly 
those that cover both Ca H&K and the Ca NIR triplet, NIR 
spectra, and good spectral sequences starting at early times 
should produce additional insight into the formation of a 
SN la SED. 

ACKNOWLEDGMENTS 

Facilities: HST(STIS) 

We thank D. Kasen, R. Kirshner, and J. Parrent for 
their comments, insights, and help. 

Based on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at the 
Space Telescope Science Institute. STScI is operated by the 



^ nr\i o r» A o A JTVTr> A o r\t\r\ FTl 



16 Foley 



R = 50 




L^^_^ R = 75 




R= 100 




R= 150 




R = 200 



3600 3700 3800 3900 
Rest Wavelength (A) 

Figure 16. Artificial Ca H&K line profiles. The profiles are the 
same as shown in Figure|3] Different resolutions (R = 50, 75, 100, 
150, and 200) are shown in each panel from top to bottom. 



Association of Universities for Research in Astronomy, Inc. 
under NASA contract NAS 5-26555. 



REFERENCES 

Altavilla G. et al., 2007, A&A, 475, 585 

Ben-Ami S., Konidaris N., Quimby R., Davis J. T., Ngeow 
C. C, Ritter A., Rudy A., 2012, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Se- 
ries, Vol. 8446, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series 

Blondin S. et al., 2006, AJ, 131, 1648 

Blondin S., Dessart L., Hillier D. J., Khokhlov A. M., 

2012a, ArXiv e-prints, 1211.5892 
Blondin S. et al., 2012b, AJ, 143, 126 

Bongard S., Baron E., Smadja C, Branch D., Hauschildt 

P. H., 2008, ApJ, 687, 456 
Branch D., Doggett J. B., Nomoto K., Thielemarm F.-K., 

1985, ApJ, 294, 619 
Branch D., Baron E., HaU N., Melakayil M., Parrent J., 

2005, PASP, 117, 545 
Branch D. et al., 2007, PASP, 119, 709 
Brown P. J. et al., 2012, ApJ, 753, 22 
Chomiuk L. et al., 2012, ApJ, 750, 164 
Chornock R., Filippenko A. V., 2008, AJ, 136, 2227 
Conley A. et al., 2011, ApJS, 192, 1 

Fabricant D., Cheimets P., Caldwell N., Geary J., 1998, 
PASP, 110, 79 

Fisher A., Branch D., Nugent P., Baron E., 1997, ApJ, 481, 
L89-h 

Foley R. J. et al., 2009, AJ, 138, 376 

Foley R. J., Brown P. J., Rest A., Challis P. J., Kirshner 

R. P., Wood-Vasey W. M., 2010, ApJ, 708, L61 
Foley R. J., Kasen D., 2011, ApJ, 729, 55 

J., Sanders N. E., Kirshner R. P., 2011, ApJ, 742, 



Foley R. 

89 
Foley R. 
Foley R. 



Foley R. J. et al. 
Foley R. J. et al. 
Foley R. J. et al., 
Garavini G. et al., 
Garavini G. et al.. 



J., 2012, ApJ, 748, 127 
J. et al., 2012a, ArXiv e-prints, 1212.2209 
2012b, ApJ, 744, 38 
2012c, ApJ, 753, L5 
2012d, ApJ, 752, 101 
2004, AJ, 128, 387 
, 2007, A&A, 471, 527 
Gerardy C. L. et al., 2004, ApJ, 607, 391 
Green J. et al., 2012, ArXiv e-prints, 1208.4012 
Hachinger S. et al., 2012, ArXiv e-prints, 1208.1267 
Hatano K., Branch D., Fisher A., Baron E., Filippenko 

A. v., 1999, ApJ, 525, 881 
Hoflich P., 1995, ApJ, 443, 89 

Hoflich P., Wheeler J. C, Thielemann F.-K., 1998, ApJ, 
495, 617 

Hogg D. W., Bovy J., Lang D., 2010, ArXiv e-prints, 

1008.4686 
Horesh A. et al., 2012, ApJ, 746, 21 
Kasen D. et al., 2003, ApJ, 593, 788 
Kasen D., Plewa T., 2007, ApJ, 662, 459 
Kelly B. C, 2007, ApJ, 665, 1489 
Kirshner R. P. et al., 1993, ApJ, 415, 589 
Lentz E. J., Baron E., Branch D., Hauschildt P. H., Nugent 

P. E., 2000, ApJ, 530, 966 



1 o O A O ATTVTTD A O f\r\f\ ITI *?•> 



SN la Line Profiles 17 

Lentz E. J., Baron E., Branch D., Hauschildt P. H., 2001, 

ApJ, 557, 266 
Li W. et al., 2001, PASP, 113, 1178 
Maguire K. et al., 2012, MNRAS, 426, 2359 
Margutti R. et al., 2012, ApJ, 751, 134 
Matheson T. et al., 2008, AJ, 135, 1598 
Matheson T. et al., 2012, ApJ, 754, 19 
Mazzali P. A. et al., 2005a, ApJ, 623, L37 
Mazzali P. A., Benetti S., Stehle M., Branch D., Deng J., 

Maeda K., Nomoto K., Hamuy M., 2005b, MNRAS, 357, 

200 

Nugent P., Phillips M., Baron E., Branch D., Hauschildt 

P., 1995, ApJ, 455, L147 
Nugent P., Baron E., Branch D., Fisher A., Hauschildt 

P. H., 1997, ApJ, 485, 812 
Nugent P. E. et al., 2011, Nature, 480, 344 
Parrent J. T. et al., 2012, ApJ, 752, L26 
Quimby R., Hoflich P., Kannappan S. J., Rykoff E., Ru- 

jopakarn W., Akerlof C. W., Gerardy C. L., Wheeler J. C, 

2006, ApJ, 636, 400 
Rodney S. A. et al., 2012, ApJ, 746, 5 
Ropke F. K. et al., 2012, ApJ, 750, L19 
Shappee B. J., Stanek K. Z., Pogge R. W., Garnavich P. M., 

2013, ApJ, 762, L5 
Silverman J. M., Ganeshalingam M., Li W., Filippenko 

A. v., 2012, MNRAS, 425, 1889 
Stanishev V. et al., 2007, A&A, 469, 645 
Suzuki N. et al., 2012, ApJ, 746, 85 
Tanaka M. et al., 2008, ApJ, 677, 448 

Tanaka M., Mazzali P. A., Stanishev V., Maurer L, Kerzen- 

dorf W. E., Nomoto K., 2011, MNRAS, 410, 1725 
Valenti S. et al., 2009, Nature, 459, 674 
Wang L. et al., 2003, ApJ, 591, 1110 
Wang X. et al, 2009a, ApJ, 699, L139 
Wang X. et al, 2009b, ApJ, 697, 380 

Williams M. J., Bureau M., Cappellari M., 2010, MNRAS, 
409, 1330 

Yaron O., Gal- Yam A., 2012, PASP, 124, 668 



18 Foley 



Table 1. Derived Quantities for lM12l Sample 



SN 




z 


Eff. Phase 




Stretch'' 


M12 vr:„>,T< 


Phase-corrected M12 


Blue JJCaHtK 


Rod DCaH&K 


Ca/Si 








(d) 








(10^ kms-')'' 


"CaH&K (10^ kl^7^)'=.= 


(10^ kms-1)'' 


(10^ kms-1)'' 


Ratio 


PTFOQdlc 


0. 


,0666 


2 


,8 


1.05 


(0.03) 


-16.99 (0.17) 


-17.78 


(0.67) 


-19.25 (0.12) 


-10.47 (0.19) 


1.709 


(0.050) 


PTF09dnl 


0. 


,019 


1 


,3 


1.05 


(0.02) 


-16.81 (0.15) 


-17.17 


(0.33) 


-18.87 (0.05) 


-9.71 (0.07) 


1.987 


(0.047) 


PTFOgfox 


0. 


,0707 


2 


,6 


0.92 


(0.04) 


-14.40 (0.04) 


-15.14 


(0.62) 


-16.99 (0.33) 


-10.28 (0.31) 


1.278 


(0.041) 


PTFOOfoz 


0. 


,05331 


2 


,8 


0.87 


(0.06) 


-14.28 (0.11) 


-15.07 


(0.67) 


-17.32 (0.30) 


-10.02 (0.28) 


1.058 


(0.006) 


PTFlObjs 


0. 


,0303 


1 


,9 


1.08 


(0.02) 


-16.47 (0.22) 


-17.02 


(0.50) 


-21.68 (0.10) 


-13.37 (0.07) 


0.790 


(0.004) 


PTFlOhdv 


0, 


,0542 


3, 


,3 


1.05 


(0.07) 


-17.03 (0.18) 


-17.94 


(0.78) 


-19.71 (0.18) 


-10.94 (0.22) 


1.637 


(0.048) 


PTFlOhmv 


0, 


,033 


2 


,5 


1.15 


(0.01) 


-15.02 (0.11) 


-15.72 


(0.59) 


-17.84 (0.16) 


-9.92 (0.18) 


1.420 


(0.026) 


PTFlOicb 


0, 


,0088 


0, 


.8 


0.99 


(0.03) 


-12.61 (0.03) 


-12.85 


(0.20) 


-16.99 (0.06) 


-10.25 (0.05) 


0.733 


(0.004) 


PTFlOmwb 


0, 


,0312 


-0. 


4 


0.94 


(0.03) 


-13.09 (0.03) 


-12.99 


(0.09) 


-16.99 (0.12) 


-10.24 (0.12) 


0.903 


(0.008) 


PTFlOqjq 


0, 


,0288 


3, 


,5 


0.96 


(0.02) 


-11.87 (0.08) 


-12.86 


(0.83) 


-17.64 (0.24) 


-10.80 (0.12) 


0.459 


(0.022) 


PTFlOtce 


0, 


,039716 


3, 


.5 


1.07 


(0.02) 


-16.02 (0.05) 


-17.00 


(0.81) 


-19.81 (0.17) 


-12.35 (0.16) 


0.850 


(0.016) 


PTFlOwnm 


0, 


,0645 


4, 


.1 


1.01 


(0.03) 


-12.90 (0.07) 


-14.04 


(0.96) 


-17.46 (0.40) 


-10.66 (0.24) 


0.579 


(0.034) 


PTFlOxyt 


0, 


0484 


3, 


.2 


1.07 


(0.04) 


-15.01 (0.08) 


-15.90 


(0.74) 


-18.74 (0.48) 


-11.62 (0.38) 


0.810 


(0.026) 


SN20091e 


0, 


,01703 


0, 


.3 


1.08 


(0.01) 


-16.08 (0.12) 


-16.16 


(0.14) 


-19.98 (0.10) 


-12.24 (0.10) 


0.966 


(0.001) 



'^'Effective phase is the measured phase divided by the stretch. 
^"As reported by M12. 

"^Measured f Ca H&K corrected by | M 1 2| velocity gradient of 280 it 230 km s""*" d 
"^Assuming a rest wavelength of 3945 A. 



Table 2. SYNOW Model Parameters for SNe 2010ae and 2011fe 



Parameter 


ol 


Nal 


Mgll 


Sill 


SlI Ca 


II 


HV Call 


CrII 


Pelll 


Coll 


SN 2011fe 


Sill A3858 


T 


0.2 


0.8 


0.5 


3 


1.2 


7 





60 


0.5 


0.4 




5 


1 


2 


2 


1.5 


2.5 




2 


2.5 


2 




8 


8 


5 


6 


12 


18 




10 


10 


10 


HV Ca 


T 


0.2 





0.5 


3 


1.2 


4 


1.7 


60 


0.5 


0.4 




5 




2 


2 


1.5 


3.5 


2 


2 


2.5 


2 


^exc 


8 




5 


15 


12 


18 


18 


10 


10 


10 


SN 2010ae 


Sili A3858 


T 


0.2 


0.8 


0.5 


1 


1.2 


4 





60 


0.5 


0.4 




5 


1 


2 


2 


1.5 


2.5 




2 


2.5 


2 




8 


8 


5 


6 


12 


18 




10 


10 


10 


HV Ca 


T 


0.2 





0.5 


1.5 


1.2 


4 


1 


60 


0.5 


0.4 




5 




2 


2 


1.5 


2 


2 


2 


2.5 


2 




8 




5 


15 


12 


18 


18 


10 


10 


10 



